The Worldsheet Formulation as an Alternative Method for Simulating 

Dynamical Fermions. 



On 



Hugo Fort 

Instituto de Ffsica, Facultad de Ciencias, 
Tristan Narvaja 1674, 11200 Montevideo, Uruguay 

The recently proposed worldsheet formulation of lattice fermions is tested for the 
first time carrying out a simulation for the simplest model: the one-flavor, strictly massless 
lattice Schwinger model. A main advantage of this alternative method for simulating 
£ — t dynamical fermions consists in its economy: it involves many fewer degrees of freedom 

. than the ordinary Kogut-Susskind formulation. The known continuum limit is reproduced 

0^ ' by the method for relatively small lattices. 

It is well known that a complete description of Quantum Chromodynamics (QCD) requires a non- 
perturbative approach. The main available non-perturbative method is lattice gauge theory. In particular, 
the Lagrangian formulation, with discretized Euclidean time, allows to apply the powerful numerical 
method of Monte Carlo simulation. However, simulating QCD with dynamical fermions or full QCD 
is still too expensive in computer time. Throuhout the years the conclusion has been the same: either 
more machine power and/or a real improvement in algorithms is needed to produce reliable estimates of 
| , hadronic quantities 

The root of the problems of QCD can be traced to the fact that the theory is denned in terms 
\ of local fields, quarks and gluons, but the physical excitations are extended composites: mesons and 

baryons. A Hamiltonian lattice formulation directly in terms of string- like excitations gave rise to the 
so-called F '-representation g]. The states |P > of this representation are described by sets of lattice 
paths associated to the extended excitations. We propose to explore a new numerical approach based 
on the recently introduced worldsheet formulation of lattice gauge theories with dynamical fermions Q] . 
This formulation, in terms of the worldsheets of the extended excitations, is the Lagrangian counterpart 
of the P-representation. The worldsheet partition function of lattice QED with staggered fermions is 
expressed as a sum over surfaces with border on self-avoiding fermionic loops. Each surface is weighted 
with a classical action written in terms of integer gauge invariant variables. The surfaces correspond 
to the worldsheets of loop-like pure electric flux excitations and meson-like configurations (open electric 
flux tubes carrying matter fields at their ends). This description, besides the general advantage of 
*^ geometric transparency, is appealing because it involves fewer degrees of freedom than the ordinary 

, Kogut-Susskind formulation implying a substantial economy of CPU time. Furthermore, the present 

formulation of dynamical fermions in terms of self-avoiding loops is closely connected with the polymer- 
like representation |^] and this might enable to exploit techniques used in other branches of physics, like 
condensed matter physics and cosmic string physics || . 

In this letter we present, for the first time, a test of the formalism developed in The simplest 
lattice gauge theory with dynamical fermions, the Schwinger model or (1+1) QED, is chosen and a Monte 
Carlo simulation is performed. This massless model can be exactly solved in the continuum. However, it 
is rich enough to share many features with 4-dimensional QCD. For this reason it has been extensively 
used as a laboratory to analyze the previous phenomena and has also become a popular benchmark to 
test different techniques to simulate theories with fermions. Our main point is to show that the method 
works properly and it might provide a novel alternative approach to simulate full QCD. 

The P-representation offers a gauge invariant description of physical states in terms of kets | P >, 
where P labels a set of connected paths P% oriented from the even sites x to the odd ones y in a lattice 
of spacing a. These paths correspond to string-like "electromeson" excitations of staggered fermions 
connected by tubes of electric flux. The internal product of a ket | > with one in terms of fields is 
given by 

$(p|) ^< py | v4 V* o„ >= ^W(py)My), (i) 

where U^{x) = exp[i0 A1 (x)] = exp[ieaA^(x)] and u and d denote respectively the up and down parts of 
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the Dirac spinor. Therefore, the lattice paths P% start in sites a; of a given parity and end in sites y with 
opposite parity. The one spinor component at each site can be described in terms of the Susskind's x(x) 
single Grassmann fields ||. The path creation operator $q in the space of kets {| P >} of a path with 
ends x and y is defined as 

<I Q = $?{x)U{Q v x )x{y). (2) 

Its adjoint operator $q acts in two possible ways annihilating the path Q v x or joining two existing 
paths in | P > one ending at x and the other starting at y. 

The worldsheet partition function Zp is a sum over the worldsheets of strings or paths of the P- 
representation. That is, surfaces Spc such that: (I) their borders T c are self-avoiding polymer-like loops 
and (II) when intersected with a time t — constant plane they produce paths beginning at even sites 
and ending at odd ones. This description at first sight is similar to the one obtained by integrating the 
fields in the Kogut-Susskind partition function Zks, nevertheless it differs in two features. In first place, 
the integration of the fields in Zks produces, besides surfaces bounded by fermionic loops, isolated links 
traversed in both opposite directions or "null" links. These "null" links were ruled out from the Zp and 
with them the myriad of different configurations corresponding to a given configuration of worldsheets. In 
second place, by virtue of the constraint (II), the surfaces of Zp when intersected with a time t = constant 
plane produce only paths beginning at even sites and ending at odd ones instead of the paths having 
ends of any parity (as it happens with the surfaces obtained from Zks)- These two differences translate 
respectively in: computer time economy and a cure for the additional species doubling problem of the 
Kogut-Susskind action. Concerning the last point, the worldsheet action has only 2 fermion species in 
4 dimensions (instead of 4 fermion species) which is phenomenologically more satisfactory. In ref. Q it 
was proved that Zp leads to the QED Hamiltonian using the transfer matrix procedure. The expression 
for Zp is as follows Q: 

z P = 5>(.nexp{-^ J2 ( 3 ) 

where n p is an integer variable attached to plaquettes (a 2-form) and a? is a sign given in terms of purely 

geometric quantities of the fermionic loops T c In D = 2, it turns out that a(J- c ) = (— l) N ^ c f ^ A ? c 
- where Nj^c, L^c and Ajra are, respectively, the number of connected parts, the length and the area of 
T c - and all the non- vanishing contributions have ajr = +1. The reason is that Npc — — jp- + Aj?c = I£, 
the number of enclosed sites by the fermionic loops T c which is always even by virtue of the above 
constraints (I) and (II) ( see ref. [f| for more details), so that we omit this factor. 

The fermionic paths T c can be expressed in terms of integer 1-forms -attached to the links- / with 
three possible values: and ± 1 with the constraint that they are non self-crossing and closed (df = 0) 
as 

^Schwingcr = £ £ exp{ __L M 2 }§if _ dn) _ (4) 

« / 

The lattice chiral condensate per-lattice-site is defined as < \X >— JFT ^2x(~^-) Xl < K^( x )>x( x )] >> 
where N s is the number of lattice sites. The corresponding operator is realized in the P-representation 
and thus we get for the chiral condensate 0] : 

<xx>= 2~lv7' (5) 

where Np is the number of connected paths at a given time t. Thus, equation (||) allows to calculate 
directly the chiral condensate simply by counting the number of "electromesons" we have when we 
intersect their world sheets with each time slice t. 
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To generate the worldsheets we use a Metropolis-type updating algorithm with the Boltzmann weight 
proportional to exp{ — ^ ^ p n^}. We simulate the model on L x L square lattices with periodic boundary 
conditions (pbc). All the plaquettes p belonging to a surface, by virtue of the self-avoiding constraint 
(I) on their borders, must fulfill the condition that the difference between the integers n p of contiguous 
plaquettes has 3 possibilities: A M n p = ±1 or (/.i = 1, 2). On the other hand, the constraint (II) implies 
that Ain p — lor — 1 according to the parity of the spatial coordinate x\. With these simple rules the 
algorithm generates the surfaces. In the computation of the chiral condensate per-lattice-site < XX > 
there is a clear difference between a simulation in the ordinary representation, in terms of fields, and 
one in the P-representation. In the first case, for the massless model, given enough time, the system 
rotates through all the degenerate minima so that < XX >= 0- Therefore, one has to calculate this order 
parameter for several small masses m/0, which select the 9 = vacuum, and then extrapolate to the 
limit m — > 0. On the other hand, in the gauge invariant P-representation from the very beginning 9 = 
JhJ and the computation with m = can be performed directly. To compute < XX > using (||) we have 
to count Np at each one of the L time slices. Therefore, for each lattice sweep we collect L values of 
< XX >• Typically, 100000 sweeps per point were performed. 

The first thing we have checked is that we get the right strong coupling behavior, both for the 
ground-state energy density ojq and the chiral condensate per-lattice-site < XX >■ The agreement with 
the series expansion is pretty good up to /3 = 0.12. 

On the other hand, the continuum theory is reached at zero lattice coupling, in the same way as 
four-dimensional asymptotically free theories like QCD. It is known that the chiral condensate in the 
continuum is given exactly by 



< ipip > 



2tt 3 / 2 



= 0.15995 (7 is Euler's constant). 



(6) 



As the lattice size L increase, the convergence of q =< xx >/ e to its known continuum value (gj) improves. 
A property of the action of (^) is that q becomes stabilized at its continuum value well inside the weak 
coupling. Before reaching this regime, this observable oscillates strongly with /3. In FIG.l we plot q vs. 
(3 for lattices of sizes ranging from L = 20 to L = 32. For L = 32 it is apparent that a value close to the 
exact continuum one is reached. 
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FIG. 1. q =< XX >/ e f° r different lattice sizes L = 20 (diamonds), L = 24 (triangles up) , L = 28 (triangles 
down) and L = 32 (circles). The horizontal line correspond to the known continuum value ~ 0.16 (statistical 
errors are indicated by vertical bars for L = 24 and L = 32). 
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Our aim has been to show that the recently introduced worldsheet formulation is a valuable alter- 
native tool in order to do numerical computations with dynamical fermions. So, at this stage, we have 
chosen the simplest algorithm: the standard Metropolis one. Obviously, there are several more sophisti- 
cated techniques that can be applied in order to improve the results. Additionally, the code can be, of 
course, optimized. 

The method presents the following advantages: 1) Easiness of computation. For instance, the lattice 
chiral condensate is diagonal in the P-representation and hence to evaluate it one simply has to count the 
number of connected open string-like excitations. Furthermore, one can compute directly < \X > f° r the 
case of strictly mass. 2) Economy I: it involves much less degrees of freedom that the Kogut-Susskind 
action. 3) Economy II: no gauge redundancy. 4) As a consequence of the constraint (II), the worldsheet 
action does not suffer from the additional species doubling problem of the Kogut-Susskind action. 

Concerning the results: we get the right continuum limit (the weak coupling fixed point) as one can 
see from the chiral condensate which is consistent with the exact results for the continuum theory. This 
is a numerical confirmation of the redundancy of the information carried out by the "null" links and the 
correctness of their elimination from the partition function. 

Finally, we would like to stress (once more) that our aim was not to present another solution to 
the Schwinger model, but to test a new general approach to tackle dynamical fermions. The results are 
promising, and this encourages us to employ more refined numerical methods in a following stage. To 
extend the method to D > 2 one has to compute the sign <j of each configuration. We are working on an 
algorithm to compute efficently this factor. 
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